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ABSTRACT 


The wake of an airfoil which oscillates in pure plunge mode is investigated in a 
water tunnel over a wide range of reduced frequency and amplitude. The main focus of 
this study 1s the comparison of the experimentally determined wake geometry with 
numerical results from a potential flow code. The wake vortices are visualized by two- 
color dye injection and velocity profiles are measured with LDV upstream and 
downstream of the airfoil. Wake signatures are examined with regard to thrust or drag 
generation. There 1s a good agreement between calculated and experimental data of the 
vortical wavelength. At high plunge velocities both approaches show a loss of wake 


symmetry and the emergence of a dual mode wake behavior. 
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I. INTRODUCTION 


The aerodynamics of oscillating airfoils is of interest both because of the role it 
plays in aircraft flutter problems and because of its application in bio-fluid-dynamics as a 
propulsive system of birds, insects and fishes. The thrust of a flapping airfoil was also used 
for propulsion of an experimantal catamaran by Schmidt [Ref. 1]. In both areas 
- aircraft aerodynamics and bio-fluid-dynamics - the interference of the wake with 
downstream located bodies or the interaction of wakes from two flapping airfroils may 
significantly improve or deteriorate flight performance or propulsive efficiency. Therefore 
the knowledge of the wake properties is of great importance. 

Freymuth [Ref. 2] explained the generation of propulsive vortical wake signatures 
of plunging and pitching airfoils and documented the time sequence of the vortex 
generation. Koochesfahani [Ref. 3] provided experimental results of the wake structure of 
a pitching airfoil and its dependence on amplitude, frequency and oscillation waveform. 
Freymuth, Gustafson and Leben [Ref. 4] compared physical visualizations and computed 
vortex dynamics of an airfoil in combined plunging and pitching motion. At the Naval 
Postgraduate School Platzer, Tuncer, Jones and Teng [Ref. 5, 6, 7] conducted a 
computational analysis of flapping airfoil aerodynamics with both a Navier-Stokes solver 
and a potential flow code. 

The objective of this investigation is to provide comprehensive experimental data 
about the wake of a flapping airfoil for the purpose of comparison with the computed 
wake structures. The wake is investigated over a wide range of frequency and amplitude. 
It is visualized in a water tunnel using dye injection technique. Velocity profiles are 


measured upstream and downstream of the airfoil using laser Doppler velocimetry. 








Il. EXPERIMENTAL APPARATUS 


A. WATER TUNNEL 


The experiments were carried out in the Naval Postgraduate School water tunnel 
facility, by Eidetics International, Inc., Model No 1520. This is a closed circuit, 
continuous flow facility with a contraction ratio of 6:1 and horizontal orientation. The test 
section 1s 38cm wide, 51cm high and 150cm long. Glass side walls, rear wall and bottom 
permit maximum optical access. There is no sealed top, thus providing simple access to 
the model. The side walls of the test section slightly diverge to compensate for boundary 
layer growth and to maintain uniform flow velocity. The flow velocity can be set in a 
range from 0 to 0.5 m/s. It is measured by a flow meter at the end of the test section and 
digitally displayed. Figure 1 depicts the layout of the water tunnel. 

For flow visualization water soluble food coloring is used. Pressurized canisters 
for 6 different colors are available. The pressurization for the dye canisters is provided by 
a small compressor and a pressure regulator is used to control the pressure level. The dye 
is injected either through dye ports in the model or with a rake upstream of the model. 
The rake provides two individually routed lines with adjustable horizontal and vertical 


position. 
B. SHAKER 
A shaker is mounted vertically on top of the test section. A Model 113 


ELEKTRO-SEIS shaker from APS Dynamics Inc. is used [Ref. 8]. This is a long stroke, 


light weight electrodynamic shaker. On its bottom end is the armature where differnt 
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Figure 1. 
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models can be attached. The shaker is operated with an APS Model 114 amplifier 
[Ref. 9] and a frequency generator which can be varied in a range from 5Hz to 60Hz. The 
amplitude can be adjusted on both the power supply (rough) and the frequency generator 
(fine). However, it cannot be set to a constant value. For a fixed amplitude setting (fixed 
current) the actual amplitude is a function of the frequency. With increasing frequency the 
amplitude decreases. This is because the unit employs a permanent magnet and is 
configured such that the armature coil remains in a uniform magnetic field over the entire 
stroke range. It 1s fundamentally a force generator. If a constant amplitude is desired, the 
force must be increased with increasing frequency. Unless the force, i.e. the current, is 
increased the amplitude decreases with increasing frequency. The magnitude of the 
generated force is directly proportional to the instantaneous value of the current. Since the 
waveform of the current is sinusoidal, the waveforms of acceleration, velocity and 


displacement are sinusoidal as well. 


C. AIRFOILS 


Two different airfoils were used: most flow visualization experiments were 
conducted with a 1cm chord steel airfoil resembling a NACA 0010. For the LDV 
measurements, a wooden NACA 0012 airfoil with a chord of 10cm was used. Both 
airfoils have a wing span of 37 cm and stretch across the whole test section. The larger 


airfoil also has dye injection ports on both surfaces. Both airfoils are shown in Figure 2. 





Figure 2. Airfoils c = 10mm and c = 100mm 


D. LASER 


The LDV measurements are performed with an OMNICHROME Model 543- 
300A Argon Ion Laser with a maximum output of 300 mW. The OMNICHROME Model 
160 power supply is used. The multiline argon beam is color-separated by a TSI Model 
9201 Colorburst Multicolor Beam Separator [Ref. 10]. The three lines generated are 
green (514.5nm), blue (488.0nm), and violet (476.5nm). Frequency shift can be added to 
each color to allow measurements in reverse flow. Thus the beam separator outputs six 
laser beams. The TSI LV Frequency Shifter Model 9186 A-4 ts used. It offers upshift and 
downshift from 2 KHz to 1OMHz. Shifted and unshifted beams are transmitted by a TSI 
Model 9275 two-component fiberoptic probe system [Ref. 11], consisting of four Model 
9271-D input couplers, four single mode polarization-preserving transmitting fibers, one 
multimode receiving fiber, the probe head and a junction box. In the probe head, the 
beams are collimated and focused with a transmitting lens. The beam separation is 50mm 
and the focal length is 350mm. The backscatter receiving optics focuses the scattered light 
onto the aperture of a TSI Model 9162 photomultiplier which is supplied by a TSI Model 
916503 photomultiplier power supply. The signal processing is done by an IFA 550 
Intelligent Flow Analyzer [Ref. 12] and the data are processed and displayed using the 
TSI Flow Information Display (FIND) software version 3.5 [Ref. 13]. The fiberoptic 
probe head is mounted on a LINTECH Model 41583 traverse table which allows 
movements in x-, y- and z-direction with an accuracy of 0.1mm. It is driven by an Applied 
Motion Products System 1618 supply unit. Figure 3 shows the water tunnel test section 
with shaker, traverse table and laser probe. 

Since the existing setup provides only one photomultiplier, only one component 


measurements (in x-direction) using the blue beams were performed. 
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Figure 3. Water Tunnel test section with shaker, 
traverse table and laser probe 





Hil. EXPERIMENTAL PROCEDURE 


A. FLOW VELOCITY 


The flow velocity in the water tunnel is set with a dial which allows settings from 
0.00 to 9.99 with an accuracy of 0.01. This setting corresponds to a velocity range from 0 
to approximately 0.5 m/s. The velocity is measured by a turbine flow meter at the end of 
the test section and digitally displayed in units ft/s. This flow meter has a poor sensitivity 
at low flow velocities. It measures the free stream velocity about 150 cm downstream of 
the position where the airfoils are mounted to the shaker. The actual flow velocity in the 
test section may be somewhat different, especially when large models or airfoils are used. 

Prior to the experiments a calibration curve was set up which relates the velocity 
setting, the displayed velocity and the actual velocity in the test section. The calibration 
data are summarized in Table 1 and the calibration curve is shown in Figure 4. The 
calibration data slightly vary with deviations of the voltage that drives the pump and with 
the amount of water filled into the tunnel. Therefore all velocity data are backed up by 


LDV measurement. 
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Table 1. Water tunnel calibration data 
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Figure 4. Water tunnel calibration curve 
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B. FREQUENCY 


The frequency at which the shaker is operated is set on an analog dial on the 
frequency generator with a resolution of 1 Hz and a range from 5 Hz to 60 Hz. All 
experiments were carried out with frquencies between 5 Hz and 10 Hz. Larger 
frequencies produce very small vortices which make it difficult to measure the wake 
geometry. Furthermore, the wake patterns tend to become blurred at high frequencies. 


All results are presented in terms of the reduced frequency k which is defined as 





where 

w@ = 2af 

f = frequency in Hz 

c =chordlength of the flapping airfoil 
U = free stream velocity in the tunnel 


High reduced frequencies are obtained with low flow velocities. 


C. AMPLITUDE 


The amplitude is set with a dial on the frequency generator which ranges from 1 to 
10. The amplitude is measured from the writing of a fixed position pencil on a paper 
which is attached to the armature of the shaker. Later, this method was improved by using 
a sharp steel needle and a wax coated paper. Thus an accuracy of + 0.25 mm = 0.025 h/c 
can be obtained. For a fixed amplitude setting, the amplitude was measured at least five 


times and a mean value is calculated. It is assumed that the airfoil performs the same 
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movement as the shaker armature. For the 1 cm chordlength airfoil this is true only for 
low frequencies (5 to 10 Hz). At higher frequencies the thin airfoil starts to bend and 
vibrate. 

For a fixed amplitude setting the actual amplitude is a function of the frequency as 
explained in section II.B. Therefore constant amplitude measurements had to be done at a 
constant frequency. 


All graphs and tables show the dimensionless amplitude h/c. 


D. DYE INJECTION 


Most experiments were done with one dye injection tube placed upstream of the 
flapping airfoil, horizontally leveled with the airfoil's zero position. The distance between 
the dye tube and the airfoil was varied in a range from a few milimeters to about 2 cm so 
that clear wake signatures were produced at the different flow velocities and amplitudes. 
Vortices originating from the lower surface during the downward motion and vortices 
originating from the upper surface during the upward motion can be clearly distinguished 
by their direction of rotation. 

Some experiments were done with two dye injection tubes with red and green dye 
in order to visualize the clockwise and the counterclockwise rotating vortices in different 
colors. This procedure was successful as well, (see Appendix A, Figure 7A through 11A 


and 18A) however, adjusting the tubes was much more difficult and time consuming. 
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E. PHOTOGRAPHY 


The vortical wake patterns are photographed with a 35 mm Nikon 2000 camera 
with manual focusing, manual aperture control and automatic shutter speed, mounted on a 
tripod. For all photographs a 100 ASA film for colour prints is used. The light is provided 


by one or more 600 watt floodlights as required. 
_#F. VORTICAL WAVELENGTH 


After every set of measurements and photographs a ruler was placed in the test 
section downstream of the airfoil location and photographed with unchanged camera 
position. The vortical wavelength is measured from the wake photographs using the 


photograph of the ruler. Figure 5 shows an example. 





Figure 5. Measurement of the vortical wavelength 
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An accuracy of 1/32 inch = 0.8 mm = 0.08 L/c is assumed. Whenever possible the mean 
value of the first three periods of fully developed vortices is calculated. If the vortical 


wavelength changes further downstream in the wake, it is not taken into account. 


G. LASER DOPPLER VELOCIMETRY 


The available LDV equipment is simple to use and does not need calibration, 


nevertheless some important points need to be considered: 


1. Alignment of the fiberoptic input couplers 


Both laser beams entering the test section must be of equal intensity. The higher 
the intensity, the better the data rate. Both - maximum intensity and equality - can be 
achieved by proper alignment of the fiberoptic input couplers. The x-, y- and z-position of 
the fiber entry must be adjusted to the beam that exits the colour separator. The TSI 


manual [Ref. 11] describes the procedure in detail. 


a Frequency shift 


To obtain a good data rate at low flow velocities and to account for negative flow 
velocities a frequency shift of 100 KHz is applied. This means that 100 KHz is added to 
each measured Doppler frequency. The frequency shifter must be set to "100 KHz down" 
and in the Optics Configuration menu of the FIND software a frequency shift of 0.1 MHz 


must be entered. 
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3. Filter setting 


According to the expected frequencies (=Doppler frequency + frequency shift) the 
filter setting for the flow analyzer must be chosen. 10 different settings are availabel in the 


Processor Setup menu of the FIND software: 


100 Hz - 1 KHz 
300 Hz - 3 KHz 
1 KHz - 10 KHz 
3 KHz - 30 KHz 
10 KHz - 100 KHz 
30 KHz - 300 KHz 
100 KHz - 1 MHz 
300 KHz - 3 MHz 
1 MHz ~ 10 MHz 
3 MHz - 15 MHz 


Selected filter: 30 - 300 KHz. 
The software also offers an automatic filter selection mode where the optimum 


filter setting is chosen according to the frequencies received. 


4, Seeding 


Seeding was not necessary. There were sufficient natural particles in the water to 


obtain a reasonable data rate. However, it was recognized that the data rate drops rapidly 


as the amount of dye in the water increases. Best results are obtained with clear water. 
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5: Sample Size and Data Rate 


The Doppler bursts occur randomly. Therefore one has to decide how many 
samples shall be collected for one data point. The sample size is entered in the Realtime 
Histogram menue. Most measurements were done with a sample size of 200 or 400. The 
data rate is the number of valid samples collected per second. It is displayed in the 
Realtime Histogram window. The data rate varied between a few Hz for low flow 
velocities in dye polluted water up to 20 or 30 Hz for higher flow velocities and clear 
water. Thus the time for collecting the data was between a few seconds and a few 
minutes. When the airfoil is flapping at 5 Hz the velocity in the wake also varies at 5 Hz 
i.e. it has many periods within the data collection time span. Therefore the LDV 


measurements produce a mean velocity. 
6. Alignment of the traverse table 
Since the Laser precisely measures the velocity component that 1s perpendicular to 


the bisector of the angle formed by the two beams the traverse table has to be exactly 


aligned so that its front face is parallel to the test section side wall. 


i] 


ie Accuracy 


The standard deviations obtained for measurements with a stationary airfoil 
(Tables 33 through 36, Appendix B) indicate that the accuracy of the LDV measurements 
of the flow velocity 1s between three and four percent. Therefore it is assumed that the 
accuracy of the quotient u/U is between six and eight percent unless the error estimate 


based on half the least significant digits results in a bigger error. 
H. MATRIX OF EXPERIMENTS 
1. Vortical Wake Pictures 
For observation and photography of vortical wake signatures experiments were 
carried out both with constant reduced frequency/variable amplitude and constant 
amplitude/variable reduced frequency. 
a. Constant reduced frequency measurements 
Ail measurements were done at a frequency of 5 Hz. The reduced 
frequency was set by choosing the flow velocity in the water tunnel. Table 2 shows the 


selected flow velocities, the resulting reduced frequency and Reynolds number, the 


amplitude range and the number of the relevant table and figure. 
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b. Constant amplitude measurements 


Most measurements were done at a frequency of 5 Hz. One series was 
measured at a frequency of 10 Hz. The reduced frequency was varied by changing the 
flow velocity in the water tunnel. One series was measured at constant flow velocity and 
with a variable frequency between 5 Hz and 12 Hz. This was very complicated because it 
required a readjustment of the amplitude after every frequency change. This procedure is 
not recommended. Table 3 shows the selected amplitude and frequency, the range of flow 
velocity, reduced frequency and Reynolds number and the number of the relevant table 
and figure. 


Table 4 shows the 119 evaluated combinations of reduced frequency and 


amplitude in one matrix. 
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velocity j|velocity jreduced j|Reynolds |amplitude |table No | figure No 
setting u frequency |number__jh/c 
[m/s) k Re 


2.10 (0.105 | 3.04 1037__—*4.0.05-0.65 | 6 | 6A, 6B 
eo 0046 | 68a | 454) 40903-0165 | 8 | BABB 
0.800.037" | 8.6 | 365 10.03-0.65 [9 9A, 9B 
0.60 40.025 412.3, | 2474 0.03-0.65 [11 LIAB 


Table 2. Measurements at constant reduced frequency 
frequency f= 5 Hz 













amplitude |frequency j|velocity |reduced j|Reynolds |table No /figure No 
h/c f u frequency | number 
Hz cm/s k Re 


Table 3. Measurements at constant amplitude 
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first two digits: Fig.No. of measurements with constant reduced 
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ts LDV Velocity Profiles 


For the first LDV measurements the same airfoil was used as for the flow 
visualization experiments. Since the airfoil is mounted to a horizontal support strut located 
downstream only upstream velocity profiles were measured. The probe volume was 
vertically traversed 1 cm upstream of the leading edge. Velocity profiles were measured 
both with and without airfoil. When the first measurements revealed a dual mode 
behaviour of the flow, a measurement of the streamwise velocity distribution above and 
below the airfoil was added in order to find an explanation for this behavior. 

For further LDV measurements a larger airfoil with 10 cm chord was used because 

® it provided better optical access for the measurement of downstream 

velocity profiles. 

= it was very stiff and thus guaranteed that the dual mode behaviour if observed 

again was not due to airfoil bending or vibration. 

= it could be adjusted to different angles of attack. Two velocity profiles were 

measured at 10° AOA. 


Table 5 shows all performed measurements and parameters. 
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IV. RESULTS AND DISCUSSION 


A. WAKE OF THE FLAPPING AIRFOIL 


The wake photographs are shown in Appendix A. Figure 6A through 11A are each 
taken at a constant reduced frequency and show the changing wake pattern as the 
amplitude h/c is increased. Appendix A, Figure 13A through 20A are each taken at a 
constant amplitude and show the changing wake pattern as the reduced frequency is 
increased. 

All Figures have some common features which are listed below. Some effects 
observed with increasing amplitude are similar to the effects of increasing reduced 


frequnecy while others are not as will be outlined. 


6 At low amplitudes a clear vortical pattern can be observed a long way 
downstream. With increasing amplitude (see Figure 9A) or increasing reduced 
frequency (see Figure 15A) the burst location moves upstream and finally only two 
or three discrete vortex pairs remain visible. 

. The vortical pattern changes with increasing amplitude/increasing reduced 
frequency. For example see Figure 8A: at low amplitudes (h/c=0.03) the wake ts 
basically symmetric. With increasing amplitude (h/c=0.08) remnant vorticity breaks 
off and finally (h/c=0.20) an asymmetric wake pattern develops. Typical 
asymmetric patterns are also found in Figure 9A(h/c=0.20), 10A(h/c=0.16), 
11 A(h/c=0.12) and 16A(k=8.6 and k=12.3) 
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‘ The vertical position of the cw and ccw vortices changes with increased 
amplitude/increased reduced frequency. For example see Figure 8A: the cw 
vortices are indicated by green dye which is ejected above the airfoil. The ccw 
vottices are indicated by red dye which is ejected below the airfoil. Initially 
(h/c=0.05) the green vortices (cw) stay above the red vortices (ccw). With 
increasing amplitude (h/c=0.08) both kinds of vortices get aligned horizontally and 
finally (h/c=0.12 through h/c=0.20) the red vortices (ccw) move above the green 
vortices (cw). The switching of the vertical position of the vortices can also be 
seen if only one color is used (for example see Figure 14A: k=3.6 cw=top, k=8.6 
cw=bottom) and is essential for the drag/thrust classification which will be 
explained later. 

° Some pictures show a downward orientation of the vortical pattern. For 
example see Figure 10A(h/c=0.03), Figure 13A(k=15.9 through k=32.2) and 
Figure 14A(k=12.3 through k=32.2). This occurs due to settling of the dye 
especially at high reduced frequencies which are obtained with very low flow 
velocities in the tunnel and at low amplitudes i.e. when there is no thrust and no jet 
velocity. 

’ The vortical wavelength i.e. the horizontal spacing between vortices with 
the same direction of rotation abviously grows as the reduced frequency becomes 
smaller (for example see Figures 13A and 18A) or as the amplitude is increased 
(see Figure 9A). 


The vortical wavelength, the drag/thrust classification and the development of an 


asymmetric wake are discussed in more detail in the following sections B through G. 
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B. VORTICAL WAVELENGTH AS FUNCTION OF AMPLITUDE 


The experimental data from Table 6 to Table 11 are plotted as function of the 
amplitude and presented in Appendix A (Figure 6B through 11B). Figure 12 shows the 
vortical wavelength as function of the amplitude for the different reduced frequencies on 
one plot. For small reduced frequencies the vortical wavelength is a linear function of the 
amplitude. For higher reduced frequencies the linear relation only holds at large 
amplitudes. As the amplitude decreases the vortical wavelength rapidly drops towards 
zero. There is a maximum amplitude at which a vortical wake pattern can be observed. If 
this amplitude is exceeded the picture turns to extreme turbulence and reverse flow as 
observed in the water tunnel. The maximum amplitude decreases with increasing reduced 


frequency as shown in Table 12. 





Table 12. Maximum amplitude for vortical wake 
formation 
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Figure 12. 


L/c as Function of Amplitude and Reduced Frequency 
Zz. 


20 


ee 


1.0 


0.5 


Oz 
Amplitude h/c 





Vortical wavelength as function of the amplitude 
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c VORTICAL WAVELENGTH AS FUNCTION OF FREQUENCY 


The experimental data from Table 13 to Table 20 are plotted as function of the 
reduced frequency and presented in Appendix A (Figure 13B through 20B). Figure 21 
shows the vortical wavelength as function of the reduced frequency for four different 
amplitudes on one plot. It shows that the vortical wavelength is inversely proportional to 
the reduced frequency. Figure 21 also shows the graph for 2x/k which represents the 
characteristic length for a periodic disturbance in a steady flow. It obviously is the limiting 
value for the vortical wavelength as the amplitude approaches zero. With increasing 
amplitude the deviation from the 27/k-line and thus the vortical wavelength are increasing 
which corresponds to the results from the measurements at constant reduced frequency. 


Similar results have been reported by Koochesfahani [Ref. 3] for a pitching airfoil. 


D. DOWNSTREAM VARIATION OF THE VORTICAL WAVELENGTH 


A careful evaluation of the wake pictures reveals that the spacing between two 
vortices of same direction of rotation is not a constant. In some cases it changes as the 
vortices move downstream. Table 21 shows a comparison between initial and final 
wavelength of selected wake pictures. The wavelength is measured from the photograph 
and the ratio of the initial and the final (downstream) value 1s calculated. It can be seen 
that as the vortices travel downstream the wavelength becomes smaller in thrust type 


wakes and it becomes larger in drag type wakes. This behaviour 1s due to viscosity and the 
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L/c as Function of Reduced Frequency and Amplitude 
8 


Amplitude 


—¥— 0.40c 
—@— 0.12c 
—a— 0.08c 
—~— 0.05c 
-~m- 2Pi/k 


O 
~—, 
» 
io 

ton) 
& 

oD) 
> 
= 
od 

O 
= 

es 

© 
> 


Reduced Frequency k 


Figure 21. Vortical Wavelength as function of the reduced frequency 
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difference between wake velocity and free stream velocity. In a thrust type wake, friction 
gradually slows the wake velocity down and the vortices are compressed. In a drag type 


wake the velocity slowly increases until it reaches free stream velocity and the vortices are 


amplitude }3 Linitial }3 Léinal 
(vortices 
No 1to4) 
mm mm 


oi | 6s foos ji4o [ios [attoia [iis [drag 
oot fos foi2 figs iso {iitor4 |o.92 [thrust 
os {a6 _|oos oso [100 [asto2e [1.25 (drag 
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Table 21 Downstream variation of the vortical wavelength 
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E. DRAG/THRUST SIGNATURE 


Thrust type and drag type of the vortex street can be clearly distinguished. The 
thrust type shows a top row of counterclockwise rotating vortices and a bottom row of 
clockwise rotating vortices, both inducing a downstream velocity component on each 
other. The drag type vortices are positioned the other way round, inducing an upstream 
velocity component on each other and eventually forming a vortex street similar to the 
Karman vortex street. Both regimes are separated by a signature where both kinds of 
vortices are horizontally aligned. The position of thrust, drag and neutral type wake 
signature as a function of reduced frequency and amplitude is shown in Figure 22. The 
separating line between the drag and the thrust regime is governed by the relationship 


k me =const 
Cc 


F. COMPARISON WITH CALCULATION 


Jones [Ref. 7] conducted computations of vortical wakes, using the unsteady 
potential flow code UPOT developed at the Naval Postraduate School by Teng [Ref. 6]. 
It is a panel code that computes and plots inviscid, incompressible flows about arbitrary 
bodies performing defined plunging or pitching motions. From the wake plots the vortical 
wavelength is measured. Thirty wake plots were analyzed at reduced frequencies 
corresponding to the experiment (k = 3.0, 4.1, 6.8, 8.6, 10.1, 12.3) and at five amplitudes 
(h/c = 0.025, 0.05, 0.1, 0.2, 0.3). The evaluated wake plots are attached in Appendix C. 
The vortical wavelengths measured from the plots are shown in Table 22 (Appendix C) 
and added to Figures 6B through 11B (Appendix A). If the calculated wavelength 1s 
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extrapolated towards h/c = 0 it intersects the axis at L/c = 2a/k. This is in agreement with 
the results from section C. A comparison between calculated and experimental values of 
the vortical wavelength as function of the amplitude and the reduced frequency is shown 
in Figure 23 and 24. There is a good agreement between computed and measured values. 
This is surprising because the experiments were done at low Reynolds numbers 
(Re=250...1000). Obviously the vortical wake can be essentially treated as inviscid flow. 

At low amplitudes and low frequencies the agreement between computed and 
measured values deteriorates and the expemmental results show a shorter vortical 
wavelength than the computed results. With decreasing amplitude or decreasing frequency 
the thrust producing wake changes into a neutral wake and finally into a drag producing 
wake, where the wavelength becomes smaller than 22/k. However, the potential flow 
code assumes inviscid flow and cannot predict drag. Therefore it predicts a larger vortical 
wavelength. 

A direct comparison between the photographs from the water tunnel experiments 
and the digital images from the panel code is shown in Figures 25 and 26 for k = 3.0 and k 
= 6.8 respectively and with h/c = 0.2, 
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Figure 24. 


L/c as Function of Reduced Frequency 
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Vortical wavelength as function of the reduced frequency, 
comparison between measurements and theory 
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Figure 25. |§ Wake comparison for k=3.00 and h/c=0.2 





Figure 26. | Wake comparison for k=6.83 and h/c=0.2 
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G. DUAL MODE BEHAVIOR 


At low frequencies and low amplitudes the wake is symmetric about a horizontal 
centerline and the wake vortices are evenly spaced. As the frequency or the amplitude is 
increased the wake symmetry is lost. The vortices begin to form pairs and travel away 
from the centerline. As the vortex pairs travel up or down, small remnants of vorticity 
break off from the vortex pairs and travel down or up. Both upward and downward travel 
are observed. Figure 27 shows both modes for k = 6.8 and h/c = 0.2. The occurrence of 
both modes is observed randomly. It is possible to make the vortices switch from mode 
"up" to mode "down" and vice versa by generating a major upstream disturbance. 
According to Jones et al [Ref. 14] the dual mode behavior is also observed in the 
computed wake plots and the direction of travel is determined by the starting condition of 
the airfoil motion. All wake photographs are evaluated with respect to a dual mode 
behavior and the result is added to Tables 6 through 11 and 13 through 20 (Appendix B) 
and Figure 22. All wakes show a dual mode behavior at high values of h-k, i.e. at high 


plunge velocities. 
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Figure 27. | Dual mode behaviour of the vortical wake 


top: mode "up", bottom: mode "down" 
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H. VELOCITY PROFILES 
1. Small Airfoil (1 cm chord) 


a. Measurement No I (upstream velocity profile) 

Figure 28 shows the velocity profile 1 ett upstream of the leading edge of 
the flapping airfoil at k = 35, h/c = 0.12. and the free stream velocity measured without 
airfoil. It shows a peak (u=1.5U) as expected for a thrust producing airfoil, however it 
also shows a minimum of 0.7 U_ which was not expected. The location of the maximum 
and the minimum are displaced which may be due to the dual mode asymmetric wake 
behavior, which is to be expected at this frequency/ amplitude Since the maximum is 


located at x =-15mm the mode 1s obviously "down". 


0 5() 


Vertical position [mm] 





Figure 28. Velocity profile 10mm upstream of leading edge 
c=10mm, h=+0.12c, f=5Hz, k=35, U =8.5mm/s 
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b. Measurement No 2 (velocity as function of the amplitude) 

Figure 29 shows the velocity upstream of the airfoil as a function of the 
amplitude, measured at the location where the velocity showed a maximum in 
measurement No | (x = -10 mm, y = -15 mm). At h/c = 0.35 a maximum velocity of 2.44 
U is attained. When the amplitude is further increased, the peak velocity drops again. The 


measurement was stopped at h/c = 0.65 because of airfoil bending vibration. 
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Figure 29. Velocity 10mm upstream of flapping airfoil leading edge 
as function of the amplitude h/c. 


c=10mm, y=-15mm, f=5Hz, k=35, U=9mm/s 
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e Measurement No 3 (velocity above and below of the airfoil) 

Figure 30 shows the streamwise velocity distribution 4mm above and 4mm 
below the stationary airfoil. It shows a small increase in velocity above (1.05 U) and 
below ( 1.1 U) the airfoil. The difference between upper and lower surface may be due to 
a small negative angle of attack. It also shows a slight minimum of 0.98 U one 
chordlength upstream of the airfoil. Because of its small magnitude this minimum is not an 


explanation for the minimum velocity observed in measurement No 1. 


Legend 


—e— below airfoil 
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Figure 30. Streamwise velocity distribution above and below of the 


stationary airfoil, c=10mm, y=+4mm, U=10.5cm/s. 
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d, Measurement No 4 (dual mode behavior) 

Like m measurement No 1 the velocity profile was measured 1 cm 
upstream of the flapping airfoil, however the amplitude was increased to h/c = 0.35 which 
is according to measurement 2 the amplitude that produces maximum flow velocities. 
Figure 31 shows the data plotted from Table 27. The measurement started at y = +100 
mm and the y-coordinate was decreased in 10 mm steps. There are two discontinuities in 
the graph: at y = 80 mm and at y = -40 mm. The wake stnucture has obviously switched 
from mode "down" to mode "up" and back again. In mode "down" the measurement was 
continued by increasing the y-coordinate again until it reached its original value y = +100 
mm. 

To verify the assumption of random mode switching the measurement was 
repeated with dye injection upstream of the airfoil. Figure 32 shows the data plotted from 
Table 28. Measurements started at y = -100 mm with positive increments. Again there is a 
discontinuity at y = -40 mm. The injected dye proved that the wake has switched from 


mode "up" to mode "down". 
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Figure 31. 
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Vertical position [mm] 


Velocity profile 10mm upstream of flapping airfoil leading edge, c=10mm, h=0.35c, f=5Hz, k=35, U=9mm/s, 


____—. = mode “up", ----- = mode “down", numbers indicate sequence of measurement 


40 


_ 


JuWIaINSeaW JO 9oUENbas 9}BoIpUl SJaquNU ‘,UMOP, OPOU = ----- ‘,dn,, spou = ——— 


‘SMW 6=() ‘SE=% ‘ZHS=} ‘OS E'O=Y “WWQ[=d ‘ospo Buripes] jlopue Suiddey jo weaysdn wwgy] syoud AON A, 


[ur] uontsod [eonia A, 


0 





ZTE SINSI,} 


45 


Me Large Airfoil (10 cm chord) 


é. Measurements No 5 and No 6 (dual mode behavior) 

Both measurements were made to confirm the observed dual mode 
upstream velocity profiles using a different airfoil. To obtain the same reduced frequency 
as in measurement No 4 the velocity was increased to 9.0 cm/s. The data from Table 29 
through Table 32 are plotted in Figure 33. Again both modes can be clearly distinguished. 


(Mode "down" was measured at 8.6 cm/s due to a small drift of the free stream velocity.) 
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Figure 33. Velocity profiles 10mm upstream of airfoil leading edge, 
c=100mm, h=0.04c, f=5Hz, k=35 (36.5 for mode "down"), 
U=9.0cm/s (8.6cm/s for mode "down") 
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f Measurements No 7 and No 8 (symmetric wake) 

Downstream velocity profiles were measured in addition to the upstream 
velocity profiles. The free stream velocity was increased to 21 cm/s to obtain a reduced 
frequency k = 15 which was small enough so that the wake no longer showed a dual mode 
behaviour. Figure 34 shows the upstream velocity profiles and Figure 35 shows the 
downstream velocity profiles both for the stationary ~ flapping airfoil. While there is 


only a small upstream effect, a jet velocity of 2U is observed in the wake. 
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Figure 34. Velocity profiles 10mm upstream of airfoil leading edge, 
c=100mm, h=0.04c, f=5Hz, k=15, U=21cm/s 
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Legend 
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Figure 35. = Velocity profiles 40mm downstream of airfoil trailing edge, 
c=100mm, h=0.04c, f5Hz, k=15, U=21cm/s 
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g. Measurements No 9 and No 10 (airfoil with AOA) 

Both measurements have the same parameters as measurements No 7 and 
No 8 except an angle of attack of 10°(pivot at 0.25c). Figure 36 shows the upstream 
velocity profiles and Figure 37 the downstream velocity profiles both for the stationary 
and flapping airfoil. The flow is directed downward and the wake is displaced towards 
negative y coordinates (similar to mode "down") and the maximum in the upstream 
velocity profiles occurs at positive y-coordinates. The maximum jet velocity is 1.9 U. It 
should be noted that the laser only measures the x-comonent of the velocity. So the total 


velocity may be larger since the velocity vector is directed downward. 
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Figure 36. Velocity profiles 10mm upstream of airfoil leading edge, 
AOA=102, c=100mm, h=0.04c, f=5Hz, k=15, U=21icm/s 
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Figure 37. Velocity profiles 40mm downstream of airfoil trailing edge, 


AOA=10°9, c=100mm, h=0.04c, f5Hz, k=15, U=21cm/s 
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V. CONCLUSIONS 


The wake of a flapping airfoil has been examined in a water tunnel by means of flow 
visualization using colour dye injection. Laser velocimetry was used for measurements of 
the flowfield upstream and downstream of the flapping airfoil. The wake geometry was 
evaluated. Different vortical patterns for thrust or drag generating wakes can be clearly 
distinguished. At high reduced frequencies and high amplitudes the thrust generating 
wakes show a dual mode behavior. The thrust vector is deflected from the horizontal 
direction either upward or downward. There is some indication that the mode of 
deflection depends on the starting condition of the airfoil motion. However modes could 
also be manipulated by upstream disturbances and random mode switching has also been 
observed. Especially the 1cm airfoil tended to random switching while the 10cm airfoil 
showed a more stable behaviour. The regions of drag generation, thrust generation and 
dual mode behaviour were determined as function of the reduced frequency k and the 
amplitude h/c. They are separated by lines (h/c) k = const. for a given airfoil. 

Qualitative and quantitative comparisons were made between the observed wake 
pictures and the wake structure computed by the unsteady potential flow code UPOT. 
There is a good agreement between measured and computed wake geometries. The dual 
mode wake structure is also found in both approaches. 

LDV measurements clearly showed the jet generated by a thrust producing flapping 
airfoil both upstream and downstream of the airfoil location. Mean streamwise velocity 
components of 2.5 times freestream velocity were observed at the upstream location 
(h/c=0.35, k=35) and of two times freestream velocity at the downstream location 
(h/c=0.04, k=15). The dual mode behaviour was also confirmed by the measured velocity 


profiles. 
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APPENDIX A, 


WAKE PHOTOGRAPHS AND VORTICAL WAVELENGTH PLOTS 


Figure 6A through 11A 


Figure 13A through 20A 
Figure 6B through 11B 


Figure 13B through 20B 


Wake photographs from constant reduced frequency 
measurements 

Wake photographs from constant amplitude measurements 
Vortical wavelength plots from constant reduced frequency 
measurements 

Vortical wavelength plots from constant amplitude 


measurements 
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Figure 6A. —_- Vortical Wake as Function of Amplitude h/c, 
f= 5Hz, u= 10.5 cm/s, k = 3.0 
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Figure 6A.(cont.) | Vortical Wake as Function of Amplitude h/c 
f= 5Hz, u= 10.5 cm/s, k = 3.0 
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Figure 7A. —- Vortical Wake as Function of Amplitude h/c 
f= 5Hz, u = 7.7cm/s, k = 4.1 
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Figure 7A.(cont.) | Vortical Wake as Function of Amplitude h/c 
f= 5Hz, u = 7.7cm/s, k = 4.1 
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Figure 8A. 


Vortical Wake as Function of Amplitude h/c 
f= 5Hz, u = 4.6cm/s, k = 6.8 
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Figure 8A.(cont.) | Vortical Wake as Function of Amplitude h/c 
f= 5Hz, u = 4.6cm/s, k = 6.8 
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Figure 9A. _—-Vortical Wake as Function of Amplitude h/c 
f= 5Hz, u = 3.7cm/s, k = 8.6 
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Figure 9A.(cont.)  Vortical Wake as Function of Amplitude h/c 
f= 5Hz, u = 3.7 cm/s, k = 8.6 


61 


+/-0.27 


+/-0.35 


+/-0.47 












. 2a PE b/c = 
191g pe bes wewieneereey 


h/co = 


a a aetnmenameanenentiond — a oe 


_h/e 





| 
a 3 A 
ne 
=o 
‘6 eit x ~ 
a ; 
“€ 
hol wees ae 
, 7 h/c al 
~~» 
ya 
—— fo 
vd 
i 
Lm me eK oad 
o. q 


Figure 10A. Vortical Wake as Function of Amplitude h/c 
f= 5Hz, u = 3.1cm/s, k = 10.1 
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Figure 10A.(cont.) | Vortical Wake as Function of Amplitude h/c 
f= 5Hz, u =3.1em/s, k= 10.1 
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Figure 11A. Vortical Wake as Function of Amplitude h/c 
f= 5Hz, u=2.5cem/s, k = 12.3 
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Figure 11A.(cont.) | Vortical Wake as Function of Amplitude h/c 
f= 5Hz, u = 2.5cem/s, k = 12.3 
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Figure 13A. Vortical Wake as Function of Reduced Frequency k 
Amplitude h = +/- 0.05c , Frequency f= 5Hz 
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Figure 13A.(cont.) | Vortical Wake as Function of Reduced Frequency k 
Amplitude h = +/- 0.05c, Frequency f = SHz 
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Figure 14A. = Vortical Wake as Function of Reduced Frequency k 
Amplitude h = +/- 0.08c, Frequency f = 5Hz 
Chordlength c = Icm 
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k=12.3 








Figure 14A.(cont.) | Vortical Wake as Function of Reduced Frequency k 
Amplitude h = +/- 0.08c, Frequency f = S5Hz 
Chord Length c = lcm 
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Figure 15A. Vortical Wake as Function of Reduced Frequency k 
Amplitude h = +/- 0.08c, U = 7.7 cm/s 
Chord Length c = lcm 
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Figure 16A. Vortical Wake as Function of Reduced Frequency k 
Amplitude h = +/- 0.12c, Frequency f = 5Hz, 
Chord Length c = lcm 
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Figure 16A.(cont.) | Vortical Wake as Function of Reduced Frequency k 
Amplitude h = +/- 0.12c, Frequency f = 5Hz, 
Chord Length c = Icm 
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k = 13.7 








Figure 17A. Vortical Wake as Function of Reduced Frequency k 
Amplitude h = +/-0.15c, Frequency f= 10 Hz, 
Chord Length c = lcm 
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Figure 17A.(cont.) 


10Hz 


‘15c 


= +/-0 


Chord Length c = lcm 
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Figure 18A.  Vortical Wake as Function of Reduced Frequency k, 
Amplitude h = +/-0.2c, Frequency f= SHz, 
Chord Length c = lem 
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Figure 18A.(cont.) | Vortical Wake at Reduced Frequency k = 1.99 
Amplitude h = +/-0.2c, Frequency f = 5Hz, 
Chord Length c = 1cm 
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Figure 19A._—- Vortical Wake as Function of Reduced Frequency k 
Amplitude h = +/-0.4c, Frequency f= SHz, 
Chord Length c = lcm 
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Figure 19A (cont.) 





Vortical Wake as Function of Reduced Frequency k, 
Amplitude h = +/-0.4c, Frequency f= 5Hz, 
Chord Lerigth c = lem 
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Figure 20A. Vortical Wake as Function of Reduced Frequency k, 
Amplitude h = +/-0.65c, Frequency f = 5Hz 
Chord Length c = 1cm 
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Figure 20A.(cont.) | Vortical Wake as Function of Reduced Frequency Kk, 
Amplitude h = +/-0.65c, Frequency f = “OH: 
Chord Length c = lcm 
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Figure 7B. 


L/c as Function of Amplitude for k=4.1 
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Vortical wavelength as function of the amplitude, reduced frequency k = 4.1 
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Figure 9B. 


L/c as Function of Amplitude for k=8.6 
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Vortical wavelength as function of the amplitude, reduced frequency k = 8.6 
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L/c as Function of Amplitude for k=12.3 


experiment 
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Figure 11B. Vortical wavelength as function of the amplitude, reduced frequency k = 12.3 
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Flapping Airfoil with Amplitude h=+/-0.08c 
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Figure 14B. Vortical wavelength as function of the reduced frequency, amplitude h/c = 0.08 
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Flapping Airfoil with Amplitude h=+/-0.15c 
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Figure 17B. Vortical wavelength as function of the reduced frequency, amplitude h/c = 0.15 
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Flapping Airfoil with Amplitude h=+/-0.4c 
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Figure 19B. Vortical wavelength as function of the reduced frequency, amplitude h/c = 0.40 
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APPENDIX B. TABLES 


Table 6 through 11 ~—-Vortical wavelength from constant reduced frequency 
measurements 
Table 13 through 20 Vortical wavelength from constant amplitude measurements 


Table 23 through 36 LDV velocity data 


YS 





amplitude amplitude vortical 
setting h/c — L_ {wavelength 






vortical 





Vortical wavelength as function of amplitude 
chordlength of flapping airfoil = 1cm = 0.39in 
displayed velocity = 0.3ft/s 
actual velocity = 0.105m/s 
frequency = 5Hz 
reduced frequency k = 3.005 
2a/k = 2.09 
Reynolds number = 1037 
wake signature: D = drag 
I = indifferent 
T = thrust 
T2= thrust, dual mode 
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amplitude vortical vortical wake signature 
h/c wavelength L | wavelength 
in. Lic 


amplitude 
setting 
oe OE a Ml 
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Table 7. Vortical wavelength as function of amplitude 

chordlength of flapping airfoil = 1cm = 0.39in 

displayed velocity = 0.2ft/s 

actual velocity = 0.077m/s 

frequency = 5Hz 

reduced frequency k = 4.058 

2a/k = 1.55 

Reynolds number = 760 

wake signature: D = drag 

= indifferent 

T = thrust 
T2= thrust, dual mode 


de | G2 | KOT NOD Fb | me 
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amplitude amplitude vortical vortical wake signature 
setting h/c ee L {wavelength 
Lic 


Table 8. 





Vortical wavelength as function of amplitude 
chordlength of flapping airfoil = lcm = 0.39in 
displayed velocity = 0. 1ft/s 
actual velocity = 0.046m/s 
frequency = SHz 
reduced frequency k = 6.828 
2a/k = 0.92 
Reynolds number = 454 
wake signature: D = drag 
I = indifferent 
T = thrust 
T2= thrust, dual mode 
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vortical vortical wake signature 
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Table 9. Vortical wavelength as function of amplitude 
chordlength of flapping airfoil = 1cm = 0.39in 
displayed velocity = 0.06ft/s (velocity setting=80) 
actual velocity = 0.037m/s 
frequency = 5Hz 
reduced frequency k = 8.564 
21/k = 0.73 
Reynolds number = 365 
wake signature: D = drag 

I =1ndifferent 
T = thrust 
T2= thrust, dual mode 
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amplitude amplitude vortical vortical wake signature 
setting h/c oa L_ | wavelength 
Lic 


Table 10. 





Vortical wavelength as function of amplitude 
chordlength of flapping airfoil = lcm = 0.39in 
displayed velocity = 0.04ft/s (velocity setting = 70) 
actual velocity = 0.031 m/s 
frequency = 5Hz 
reduced frequency k = 10.105 
2a/k = 0.62 
Reynolds number = 306 
wake signature: D = drag 

I =1:indifferent 

T = thrust 

T2= thrust, dual mode 
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amplitude amplitude vortical vortical wake signature 
setting h/c a L | wavelength 





Table 11. —- Vortical wavelength as function of amplitude 
chordlength of flapping airfoil = 1cm = 0.391in 
displayed velocity = 0.02ft/s 
actual velocity = 0.025m/s 
frequency = 5Hz 
reduced frequency k = 12.323 
27a/k = 0.51 
Reynolds number = 247 
wake signature: D = drag 

= indifferent 
T = thrust 
T2= thrust, dual mode 
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Table 13.  Vortical wavelength as function of reduced frequency 
chordlength of flapping airfoil = 1cm = 0.39in 
amplitude = 0.5mm = 0.05c 
frequency = 5Hz 
Reynolds number = 99...1037 
wake signature: D = drag 

I =indifferent 
T = thrust 
T2= thrust, dual mode 
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Table 14. _—- Vortical wavelength as function of reduced frequency 
chordlength of flapping airfoil = 1cm = 0.391n 
amplitude = 0.8mm = 0.08c 
frequency = 5Hz 
Reynolds number = 49...869 
wake signature: D = drag 

I = indifferent 
T = thrust 
T2= thrust, dual mode 





F —? Reduced | Vortical Vortical Wake 
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Table 15. —_ Vortical wavelength as function of reduced frequency 
chordlength of flapping airfoil = 1cm = 0.39in 
amplitude = 0.8mm = 0.08c 


velocity setting 1.66 
displayed 0.20ft/s 
actual 0.077m/s 


Reynolds number = 760 
wake signature: D = drag 
I =indifferent 
T = thrust 
T2= thrust, dual mode 
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Table 16. _—- Vortical wavelength as function of reduced frequency 
chordlength of flapping airfoil = 1cm = 0.39in 
amplitude = 1.2mm = 0.12c 
frequency = 5Hz 
Reynolds number = 138...869 
wake signature: D=drag 

I = indifferent 
T = thrust 
T2= thrust, dual mode 
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Table 17. —- Vortical wavelength as function of reduced frequency 
chordlength of flapping airfoil = 1cm = 0.39in 
amplitude = 1.5mm = 0.15c 
frequency = 10Hz 
Reynolds number = 454...2844 
wake signature: D = drag 
I = indifferent 
T = thrust 
T2= thrust, dual mode 
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Table 18. Vortical wavelength as function of reduced frequency 
chordlength of flapping airfoil = lcm = 0.39in 
amplitude = 2mm = 0.2c 
frequency = 5Hz 
Reynolds number = 454...1560 
wake signature: D = drag 

= indifferent 
T = thrust 
T2= thrust, dual mode 
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Table 19. —_ Vortical wavelength as function of reduced frequency 
chordlength of flapping airfoil = 1cm = 0.39in 
amplitude = 4mm = 0.4c 
frequency = 5Hz 
Reynolds number = 454...2844 
wake signature: D = drag 

I =aindifferent 
T = thrust 
T2= thrust, dual mode 
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Table 20. —- Vortical wavelength as function of reduced frequency 
chordlength of flapping airfoil = 1cm = 0.39in 
amplitude = 6.5mm = 0.65¢ 
frequency = 5Hz 
Reynolds number = 1560...2844 
wake signature: D = drag 

I = indifferent 
T = thrust 
T2= thrust, dual mode 
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Table 23. Velocity profile 10mm upstream of flapping airfoil leading edge, 
c=10mm, h=+0.12c, f=5Hz, k=35, sample size = 200 
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y-position [mm] _ | velocity [m/s] stand. dev. [m/s] | turb. intensity [%] 
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Table 24. Water tunnel free stream velocity profile 
' (same x-coordinate as in Table 23, without airfoil) 
sample size = 200 
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Table 25. Velocity 10mm upstream of flapping airfoil leading edge as function 
of the amplitude h/c, c=10mm, y=-15mm, f= 5Hz, U=9mm/s, k=35, 
sample size = 200. 
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Table 26. Streamwise velocity distribution above and below of the 


Stationary airfoil, c=10mm, y=+4mm, U=105mm/s 
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Table 27. Velocity profile 10mm upstream of flapping airfoil leading edge, 
c=10mm, h=+0.35c, f5Hz, k=35, U=9mm/s 
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Velocity profile 10mm upstream of flapping airfoil leading edge, 
c=10mm, h=+0.35c, f=5Hz, k=35, U=9mm/s 
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Table 29. Velocity Profile 10mm upstream of flapping airfoil leading edge, 
mode "up", c=100mm, h=+0.04c, f=5Hz, k=35, U=9cm/s, 
sample size = 400 
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Table 30. Velocity profile 10mm upstream of stationary airfoil leading edge, 
c=100mm, U=9cm/s, sample size = 200 
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Table 31. Velocity profile 10mm upstream of flapping airfoil leading edge, 
mode "down", c=100mm, h=+0.04c, f=5Hz, k=36.5, U=8.6cm/s, 
sample size = 400 
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Table 32. Velocity profile 10mm upstream of stationary airfoil leading edge, 
c=100mm, U=8.6cm/s, sample size = 400 
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Table 33. Velocity profiles 10mm upstream of airfoil leading edge, c=100mm, 
h=+0.04c, f=5Hz, k=15, U=21cm/s, sample size = 400 
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Table 34. Velocity profiles 40mm downstream of airfoil trailing edge, c=100mm, 
h=+0.04c, f=5Hz, k=15, U=21cm/s, sample size = 400 
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Table 35. Velocity profiles 10mm upstream of airfoil leading edge, AoA=10°, 
c=100mm, h=+0.04c, f=5Hz, k=15, U=21¢cm/s, sample size = 400 
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Table 36. Velocity profiles 40mm downstream of airfoil trailing edge, AoA=10°, 
c=100mm, h=+0.04c, f=5Hz, k=15, U=21cm/s, sample size = 400 
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APPENDIX C. COMPUTED WAKES 


Figure 38 trough 43 Wake plots from Ref. [7] 


Table 22 Vortical wavelength measured from the wake plots 
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April 28, 1995 


UPOTG wake geometry; & = 3.00, 60 time-steps/cycle, NACA 0010 airfoil defined by 100 points. 
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Figure 38. Wake geometry for k=3.00 
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April 28, 1995 


teps/cycle, NACA 0010 airfoil defined by 100 points. 
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Wake geometry for k=4.06 
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UPOTG wake geometry; & = 6.83, 60 time-steps/cycle, NACA 6010 sirfoi! defined by 100 points. 
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Figure 40. Wake geometry for k=6.83 
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April 28, 1995 
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UPOTG wake geometry; & = 8.56, 60 litce-sleps/cycle, NACA GGi6 atrfotl defined vy 160 potuts. 
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Figure 41. Wake geometry for k=8.56 
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April 28, 1995 


UPOTG wake geometry; & = 10.11, 60 time-steps/cycle, NACA 0010 airfoil defined by 100 points. 
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Figure 42. Wake geometry for k=10.11 
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April 28, 1995 


UPOTG wake geometry; k = emis 60 eae heats NACA 0010 airfoil defined by 100 points. 
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Figure 43. 
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Table 22. Vortical Wavelength L/c as Function of Amplitude h/c and Reduced 


Frequency k, 
measured from Panel Code computed wake plots 
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